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1. Introduction 

Experiments designed to measure the redshifted 21 cm line from the Cosmic Dawn (CD) and 
Epoch of Reionization (EoR) are challenged by the strong astrophysical foreground contamination, 
ionospheric distortions, radio frequency interference and complex instrumental response. In order 
to reliably detect the cosmological signal from the observed data, it is essential to understand in 
detail all aspects of the experiment. For example, the cosmological signal has some characteristics 
which differentiates it from the foregrounds and noise. By use of proper statistics, it is possible to 
remove these components to extract signatures of reionization (see Chapman et al. 2014). 

To alleviate some of the problems associated with the observations of the weak cosmological 
signal, several cross-correlation analyses with observations in other frequency windows have been 
proposed. In doing this, the noise/systematics in two observations of different frequencies and 
strategies may cancel out. Even if SKA may have a high enough sensitivity not to need cross¬ 
correlation techniques in order to detect the signal, cross-correlating with other probes will improve 
our understanding of the process of reionization (for an overview Mellema et al. 2013). 

In this paper we present potential of cross-correlation studies between the SKA cosmological 
21 cm data, which reflects the neutral hydrogen content of the Universe as a function of redshift, 
with: (i) the kinetic Sunyaev-ZeTdovich (kSZ) effect in the Cosmic Microwave Background (CMB) 
data, produced by the scattering of CMB photons off free electrons produced during the reionization 
process (Sec. 2); (ii) galaxy surveys (Sec. 3); and (iii) the near infrared (NIR) backgrounds that 
reflect the primordial star formation (Sec. 4). 

Throughout this study we consider three SKA configurations: early science SKA1-LOW, 
SKA1-LOW and SKA2-LOW. The SKA1-LOW (Dewdney et al. 2013) is our standard SKA con¬ 
figuration for the noise simulations, using the OSKAR simulator 1 . We assume full correlation 
between all 866 core stations of SKA1-LOW, where the maximum baseline length is 5.29 km. 
Baseline coordinates are generated for a 12-hour synthesis observation, with a 5-minute sampling 
interval. The noise is then re-normalised to reflect 1000 hours of integration. For early science 
SKA1-LOW we assume that the sensitivity is halved. For SKA2-LOW we assume that the sensi¬ 
tivity is quadrupled. 

2. The Cross-Correlation with the kSZ 

One of the leading sources of secondary anisotropies in the CMB is due to the scattering of 
CMB photons off free electrons (Zeldovich & Sunyaev 1969). The effect of anisotropies when 
induced by thermal motions of free electrons is called the thermal Sunyaev-ZeTdovich effect (tSZ) 
and when due to bulk motion of free electrons is called the kinetic Sunyaev-ZeTdovich effect (kSZ). 
The latter is far more dominant during reionization (for a review of secondary CMB anisotropies 
see, e.g. Aghanim et al. 2008). 

The kSZ effect from a homogeneously ionized medium, i.e., with ionized fraction only as a 
function of redshift, has been studied both analytically and numerically by a number of authors; 
the linear regime of this effect was first calculated by Sunyaev & Zeldovich (1970) and subse¬ 
quently revisited by Ostriker & Vishniac (1986) and Vishniac (1987) - hence also referred to as 

1 http://www.oerc.ox.ac.uk/ ska/oskar 
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Figure 1: An example of the cross-power spectrum of the kSZ and the cosmological 21cm signal at z = 11. 
The solid line is for a ‘patchy’ reionization history, while the dashed line is for a ‘homogenous’ history. 
(Tashiro et al. 2011) 


the Ostriker-Vishniac (OV) effect. In recent years various groups have calculated this effect in its 
non-linear regime using semi-analytical models and numerical simulations (Gnedin & Jaffe 2001; 
Santos et al. 2003; Zhang et al. 2004). These studies show that the contributions from non-linear 
effects are only important at small angular scales (/ > 1000), while the OV effect dominates at 
larger angular scales. 

The kSZ effect from patchy reionization was first estimated using simplified semi-analytical 
models by Santos et al. (2003), who concluded that it dominates over that of a homogeneously 
ionized medium. More detailed modeling of the effect of patchy reionization were subsequently 
performed using numerical simulations (Salvaterra et al. 2005; Iliev et al. 2007) and semi-analytical 
models (McQuinn et al. 2005; Zahn et al. 2005; Mesinger et al. 2012). Dore et al. (2007) used 
numerical simulations to derive the expected CMB polarization signals due to EoR patchiness. The 
CMB bolometric arrays Atacama Cosmology Telescope (ACT, Fowler et al. 2010) and South Pole 
Telescope (SPT, Shirokoff et al. 2011) are currently being used to measure the CMB anisotropies 
at the scales relevant to reionization (3000 < i < 8000). The SPT results are starting to put limits 
on the duration of reionization (Zahn et al. 2012). 

Cross-correlation between the cosmological 21cm signal, as measured with SKA, and the 
secondary CMB anisotropies provide a potentially useful statistic. The cross-correlation has the 
advantage that the measured statistic is less sensitive to contaminants such as the foregrounds, 
systematics and noise in comparison to “auto-correlation” studies. 

Analytical cross-correlation studies between the CMB temperature anisotropies and the EoR 
signal on large scales (/ ~ 100) were canned out by Alvarez et al. (2006); Adshead & Furlanetto 
(2008); Lee (2009) and on small scales (/ > 1000) by Cooray (2004); Salvatena et al. (2005); Slosar 
et al. (2007); Tashiro et al. (2008, 2010, 2011). Cross-correlation between the E- and B-modes of 
CMB polarization with the redshifted 21cm signal was done by Tashiro et al. (2008); Dvorkin et al. 
(2009). Numerical studies of the cross-conelation were canied out by Salvatena et al. (2005); Jelic 
et al. (2010). 

These studies showed that the kSZ and the redshifted 21cm signal: (i) anti-correlate on the 
scales conesponding to the typical size of ionized bubbles; and (ii) correlate on the larger scales, 
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Figure 2: S/N of the 21cm cross-correlation with the kSZ for two different reionization models, defined 
with the redshift at which the ionized fraction equals 0.5 (z re ) and the reionization duration (Az). The S/N 
is given as a function of the sky fraction and the normalised noise power spectrum. We assume lOOOh of 
integration time with early SKA1-LOW, SKA1-LOW and SKA2-LOW and a field of view of 25 deg 2 . For 
the CMB data we assume Planck sensitivity, (based on Tashiro et al. 2010) 

where the patchiness of the ionization bubbles are averaged out (see Fig. 1). The significance of 
the anti-correlation signal depends on the reionization scenario (Salvaterra et al. 2005; Jelic et al. 
2010; Tashiro et al. 2011). 

The cross-correlation signal turns out to be difficult to detect (see Fig. 2). We might be able to 
make detection only with SKA2-LOW, assuming a very radical reionization duration of Az = 0.01. 
If duration of reionization is longer, detection is not possible for all three SKA configurations. 
However, the kSZ signal induced during the EoR could possibly be detected in the power spectra 
of the CMB and used to place some additional constraints on this epoch in the history of our 
Universe. 

3. The Cross-Correlation with Galaxy Surveys 

Following Lidz et al. (2009), one can define the cross power spectrum between the 21cm 
emission and the galaxies as: 


A 2 


Fit,gal W — ^21,gal(^)/^^0 


*HI 




(3.1) 


where 8T^q is the 21cm brightness temperature relative to the CMB for neutral gas at the mean 
density of the universe, xm is the neutral fraction and A 2 h (k) is the dimensionless cross power 
spectrum between fields a and b. In order to construct the cross power spectrum, one therefore 
requires three fields, the density field (p), the neutral hydrogen field (xhi), and the galaxy field 
(gal), which can be obtained via numerical simulations of galaxy formation and the reionization 
process. 

It is found that the 21cm emission is initially correlated with galaxies on large scales, anti¬ 
correlated on intermediate scales, and uncorrelated on small scales. This picture quickly changes 
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Figure 3: The circularly averaged, unnormalized 2D 21cm - galaxy cross power spectrum (A 2 j gal (k); and 
correlation coefficient for two redshifts. In the analysis we followed Wiersma et al. (2013). 


as reionization proceeds and the two fields become anti-correlated on large scales (Lidz et al. 2009; 
Wiersma et al. 2013). These (anti-) correlations can be a powerful tool in indicating the topology 
of reionization and should form important diagnostic tools for SKA observations. 

If the effect of observing and selecting real galaxies is taken into account, the result depends on 
the observational campaign considered. For example, for a drop-out technique (as in observations 
of Lyman Break Galaxies), the normalization of the cross power spectrum seems to be the most 
powerful tool for probing reionization. In particular, it is quite sensitive to the ionized fraction 
as different reionization histories yield similar cross power spectra for a fixed ionized fraction. 
When instead a more precise measurement of the galaxy redshifts is available (as in Ly-a Emitters 
surveys), which provides the three-dimensional position of the galaxy, much more information 
about the nature of reionization can be extracted, since the shape and the normalization of the 
cross power spectrum provide useful information. In addition, the observability of the Ly-a line 
from these galaxies is affected by neutral patches in the IGM and thus Ly-a Emitters surveys are 
particularly useful for EoR studies (McQuinn et al. 2007; Jensen et al. 2014). 

Figure 3 shows the 21cm - Ly-a Emitters cross power spectrum (with and without SKA noise) 
and correlation coefficient for two redshifts. Here the noise assigned to the 21cm survey is the one 
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of the SKA telescope with SKA 1-LOW configuration and noise for 1000 hours of integration time, 
while the Ly-a Emitters survey has the same characteristics of the one described in Jensen et al. 
(2014); Ouchi et al. (2010) with the Subaru telescope. As we can see, the 21cm - Ly-a Emitters 
cross power spectrum is not very sensitive to noise. Therefore, early SKA 1-LOW and SKA2-LOW 
will show very similar results. The effect of neutral patches in the IGM on the observability of 
these Ly-a Emitters is not included here. 

4. The Cross-Correlation with Near Infrared Backgrounds 

Understanding star formation at high redshifts is fundamentally linked to our understanding of 
reionization history. In order for reionization to occur, a plentiful source of ionizing photons was 
needed. Stars are a likely candidate to produce these photons. Therefore, any attempt to understand 
reionization must be paired with an understanding of primordial star formation. 

However, current observational constraints suggest that in order for the Universe to be reion¬ 
ized, a plentiful number of small galaxies beneath the detection limit of current surveys arc most 
likely needed. Because they arc so faint, it is difficult to obtain information about these galaxies 
directly. Instead, these galaxies can be observed using indirect means. For example, the cumulative 
light from these galaxies would be redshifted to the infrared, and therefore, any background in the 
infrared may provide clues to the nature of these galaxies. Many have suggested that the remnant 
light from these galaxies is indeed present in the excess emission in the Near Infrared Background 
(NIRB), and if this is true, looking at the mean intensity and the power spectrum of the NIRB can 
give information about these early stellar - populations (e.g., Santos et al. 2002; Magliocchetti et al. 
2003; Salvaterra & Ferrara 2003; Cooray 2004; Kashlinsky 2005; Fernandez & Komatsu 2006; Fer¬ 
nandez et al. 2010; Cooray et al. 2012; Fernandez et al. 2012; Kashlinsky et al. 2012; Fernandez & 
Zaroubi 2013; Yue et al. 2013). 

On the other hand, 21cm emission originates from areas of neutral hydrogen, so it corresponds 
to regions that have not seen plentiful star formation and hence have not yet been ionized. There¬ 
fore, we should expect that regions that are bright in 21cm background emission should not be 
bright in the infrared. Conversely, regions that are bright in the infrared have plenty of star forma¬ 
tion and therefore should be dim in the 21cm background. Because of this, these two observations 
should be anti-correlated. 

In Fernandez et al. (2014), we predict the observational properties of both high redshift galax¬ 
ies (which would be observable in the infrared) and the neutral regions that have yet to be ionized 
by galaxies (observable as part of the 21cm background). In order to do this, we combined large 
scale N-body simulations (Iliev et al. 2014), which included radiative transfer, with analytical mod¬ 
els of the luminosities of the halos (Fernandez & Komatsu 2006). The brightness temperature of 
the 21cm background was also computed from the simulation, which depends on the density and 
neutral fraction of the IGM around the galaxies. 

In Fernandez et al. (2014), we then produced simulated sky maps for both the 21cm emission 
and the emission in the infrared. These maps were created by combining the emission models and 
the simulation output to create luminosity cubes. These cubes were then randomly rotated and 
stacked along the line of sight to create a three dimensional cuboid. This was projected onto two 
dimensions to create a simulated sky map in both the radio and the infrared. In Fernandez et al. 
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Figure 4: Left: The cross correlation of the entire NIRB against the 21cm background, assuming three SKA 
configurations with lOOOh of integration time. The results do not depend heavily on noise levels, indicating 
that errors from unknowns in observation dominate over errors introduced from the noise. Right: The cross 
correlation at various redshifts: 6.4 < z < 7, 7 < z < 8, and 8 < z < 11, assuming SKA1-LOW configuration. 


(2014), we created sky maps specific for the LOFAR instrument; here, we present maps tailored 
for three SKA configurations: early science SKA1-LOW, SKA1-LOW and SKA2-LOW (assuming 
lOOOh of integration). 

For the radio map, since the 21cm background is line emission, each frequency corresponds to 
a different redshift, and observations can be adjusted to correspond to any arbitrary redshift range. 
On the other hand, since the emission in the infrared has a considerable amount of continuum, it is 
very difficult to extract redshift information for observations at any given wavelength. Therefore, 
the map in the infrared is created using the integrated light from a redshift of 6 to 30, or the entire 
simulation volume. 

These two maps in the infrared and in the radio can then be cross-correlated using the Pearson 
correlation coefficient: 


cov((8T b ),I N i RB ) 

P21cm,NIRB — - 

GST^Inirb 


(4.1) 


where Inirb is the intensity in the infrared and 8T\ j, is the brightness temperature of the 21cm line. 
The results of this cross-correlation is shown in left panel of Fig. 4, which illustrates the cross¬ 
correlation between the integrated light from 6 < z < 30 in the near-infrared against the entire 21cm 
map from 6.4 < z < 11, assuming three SKA configurations. Error bars are generated from cross- 
correlating the 21cm emission with 1000 randomized maps in the infrared. The cross-correlation 
coefficient shows that these two maps are strongly anti-correlated, corresponding to the fact that 
the two emission maps are sensitive to emission from mutually exclusive areas. In addition, there 
is not a significant difference in result between Early SKA1-LOW, SKA1-LOW, and SKA2-LOW, 
indicating that the errors introduced from the noise is sub-dominant. 

In addition, since the 21cm background is line emission, we can cross-correlate only specific 
slices of the 21cm background with the entire NIRB. Results of this for SKA1-LOW configuration 
are shown in the right panel of Fig. 4. The correlation coefficient is computed using radio maps for 
early (z > 8), mid (6.9 < z < 8) and late (z < 6.9) stages of reionization against the entire NIRB. 
The correlation coefficient is the most negative during mid stages of reionization. At early stages 
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of reionization, ionized bubbles around the host galaxies arc too small, perhaps smaller than the 
smoothing length, and therefore no correlation is seen between the two maps. At late stages of 
reionization, most of the structure in the 21cm maps has disappeared. 

5. Conclusions 

In the view of presented results, cross-correlations with other probes will improve our under¬ 
standing of the process of reionization. The cross-correlation with the galaxy surveys and NIR 
backgrounds will specifically help in answering the question which types of galaxies are mostly 
responsible for reionization. There is not a significant difference in results between Early SKA1- 
LOW, SKA1-LOW, and SKA2-LOW, indicating that the errors introduced from the noise is sub¬ 
dominant. Unfortunately, the cross-correlation signal in the case of the CMB turns out to be difficult 
to detect. There is only a possibility of detection in the case of a very short reionization scenario 
and using SKA2-LOW. 
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